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ABSTRACT 

We present a study on the effect of undetected stellar companions on the derived planetary radii for 
the Kepler Objects of Interest (KOIs). The current production of the KOI list assumes that the each 
KOI is a single star. Not accounting for stellar multiplicity statistically biases the planets towards 
smaller radii. The bias towards smaller radii depends on the properties of the companion stars and 
whether the planets orbit the primary or the companion stars. Defining a planetary radius correction 
factor Xji, we find that if the KOIs are assumed to be single, then, on average, the planetary radii may 
be underestimated by a factor of {Xb) ~ 1.5. If typical radial velocity and high resolution imaging 
observations are performed and no companions are detected, this factor reduces to {Xr) rs 1.2. The 
correction factor {Xr) is dependent upon the primary star properties and ranges from (A/j) Ri 1.6 for 
A and F stars to {Xr) k, 1.2 for K and M stars. For missions like K2 and TESS where the stars may 
be closer than the stars in the Kepler target sample, observational vetting (primary imaging) reduces 
the radius correction factor to {Xr) r; 1.1. Finally, we show that if the stellar multiplicity rates are 
not accounted for correctly, occurrence rate calculations for Earth-sized planets may overestimate the 
frequency of small planets by as much as 15 — 20%. 

Subject headings: (stars:) planetary systems, (stars:) binaries: general 


1. INTRODUCTION 


The Kepler Mission [Borucki et al.l (120101) . with the dis¬ 
covery of over 4100 planetary candidates in 3200 sys¬ 
tems, has spawned a revolution in our understanding of 
planet occurrence rates around stars of all types. One 
of Kepler’s profound discoveries is that small planets 
(Rn 3Rm) are nearly ubiquitous (e.g.. [Howard et ahl 

nneaul[2013Hl 


201^ Dressing fc Charbonneaijl201^lFressin et al.ll2013t 

Petigura et al.ll2013l : iBatalhal 12014^ and, in particular, 

some of the most common planets have sizes between 
Earth-sized and Neptune-sized - a planet type not found 
in our own solar system. Indeed, it is within this group 
of super-Earths to mini-Neptunes that there is a tran¬ 
sition from “rocky” planets to “non-rocky planets”; the 
transition is near a planet radius of 1.6i?0 and is very 
shar p - occurring within Ri 0.2Rm of this transition ra¬ 
dius (IMarcv et al.ll20l4 IRogersll2014D . 

Unless an intra-system comparison of planetary radii 
is performe d where only the r elative planetary sizes are 
important (iCiardi et al.l 120131) . having accurate (as well 
as precise) planetary radii is crucial to our comprehen¬ 
sion of the distribution of planetary structures. In par¬ 
ticular, understanding the radii of the planets to within 
^ 20% is necessary if we are to understand the relative 
occurrence rates of “rocky” to “non-rocky” planets, and 
the relationship between radius, mass, and bulk density.. 
While there has been a systematic follow-up observa¬ 
tion program to obtain spectroscopy and high resolution 
imaging, only approximately half of the Kepler candi¬ 
date stars have been observed (mostly as a result of the 
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brightness distribution of the candidate stars). Those 
stars that have been observed have been done mostly to 
eliminate false positives, to determine the stellar param¬ 
eters of host stars, and to search for nearby stars that 
may be blended in the Kepler photometric apertures. 

Stars that are identified as possible binary or triple 
stars are noted on the Kepler Community Follow-Up 
Observation Program we bsitjj, and are often handled 
in ind ividual papers fe.g.. lStar et al.l[^0l4 lEverett et al.1 
I20I4D . The false positive assessment of an KOI (or all 
of the KOIs) can take into account the likelihood of stel¬ 
lar c ompanions fe.g.. [Morton fc JohnsonI 120111 : iMortonl 
I2012D . and a false positive probability will likely be in¬ 
cluded in future KOI lists. But presently, the current 
production of the planetary candidate KOI list and the 
associated parameters are derived assuming that all of 
the KOI host stars are single. That is, the Kepler 
pipeline treats each Kepler candidate host star as a 
single star (e.g., iBatalha et al.l 120111 : IBurke et al.ll20l4 
IMullallv et al.ri2014[) . Thus, statistical studies based 
upon the Kepler candidate lists are also assuming that 
all the stars in the sample set are single stars. 

The exact fraction of multiple stars in the Kepler can¬ 
didate list is not yet determined, but it is certainly 
not zero. Recent work suggests that a non-negligible 
fraction (^ 30 — 40%) of the Keple r host stars may 
be inultiple stars (lAdams et al.ll2012L |20l5 ILaw et all 
12014 [Dressing et al.112014 iHorch et al.ll2014D , although 
other work may indicate that (giant) plan et formation 
may b e suppressed in multiple star systems (iWang et al.1 
I2014allbl) . The presence of a stellar companion does not 
necessarily invalidate a planetary candidate, but it does 
change the observed transit depths and, as a result, the 
planetary radii. Thus, assuming all of the stars in the Ke- 

^ https://cfop.ipac.caltech.edu/ 
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pier candidate list are single can introduce a systematic 
uncertainty into the planetary radii and occurrence rate 
distributions. This has already been discussed for the oc¬ 
currence rate of hot Jupiters in the Kepler sample where 
it was found that ^ 13% of hot Jupiters were classi¬ 
fied as smaller planets because of the unacco unted effects 
of tran sit dilution from stellar companions (iWang et al.l 

In this paper, we explore the effects of undetected stel¬ 
lar gravitationally bound companions on the observed 
transit depths and the resulting derived planetary radii 
for the entire Kepler candidate sample. We do not 
consider the dilution effects of line-of-sight background 
stars, rather only potential bound companions, as com- 
panio ns within 1" are most likely bound c ompanions 
fe.g.. lHorch et al.ll20l4 IcTilliland et al.ll2ni4[ l. and most 
stars beyond 1" ar e either in the Kepler Input Catalog 
(jBrown et al.ll20nll or in the UKIRT survey of the Ke¬ 
pler field and, thus, are already accounted for with re¬ 
gards to flux dilution in the Kepler project transit fitting 
pipeline. Within I", the density of blended background 
stars is fairly low, rang ing between 0.001 — 0.007 stars/D" 
(|Lillo-Box et al.l 120141 . Thus, within a radius of I", we 
expect to find a blended background (line-of-sight) star 
only 0.3—2% of the time. Therefore, the primary contam¬ 
inant within I" of the host stars are bound companions. 

We present here probabilistic uncertainties of the plan¬ 
etary radii based upon expected stellar multiplicity rates 
and stellar companion sizes. We show that, in the ab¬ 
sence of any spectroscopic or high resolution imaging 
observations to vet companions, the observed planetary 
radii will be systematically too small. However, if a can¬ 
didate host star is observed with high resolution imaging 
(HRI) or with radial velocity (RV) spectroscopy to screen 
the star for companions, the underestimate of the true 
planet radius is signihcantly reduced. While imaging and 
radial velocity vetting is effective for the Kepler candi¬ 
date host stars, it will be even more effective for the K2 
and TESS candidates which will be, on average, 10 times 
closer than the Kepler candidate host stars. 

2. EFFECTS OF COMPANIONS ON PLANET RADII 

The planetary radii are not directly observed; rather, 
the transit depth is the observable which is then related 
to the planet size. The observed depth (Jo) of a plane¬ 
tary transit is defined as the fractional difference in the 
measured out-of-transit flux {Ftotai) and the measured 
in-transit flux {Ftransit)'- 

j. Fj^Qidi — ^transit /-, \ 

^o = - - -. (1) 

^total 

If there are N stars within a system, then the total out- 
of-transit flux in the system is given by 

N 

Ftotai ( 2 ) 

and if the planet transits the star in the system, then 
the in-transit flux can be defined as 

Ftransit = Ftotai ~ Ft {Rp/Rt*) (3) 

where Ft is the flux of the star with the transiting planet, 
Rp is the radius of the planet, and is the radius of 


the star being transited. Substituting into equation 
the generalized transit depth equation (in the absence of 
limb darkening or star spots) becomes 


<5 
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(—) 

\ Ftotai / 


RtJ 
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(4) 


For a single star, Ftotai = Ft and the transit depth 
expression simplihes to just the square of the size ratio 
between the planet and the star. However, for a multiple 
star system, the relationship between the observed tran¬ 
sit depth and the true planetary radius depends upon the 
brightness ratio of the transited star to the total bright¬ 
ness of the system and on the stellar radius which changes 
depending on which star the planet is transiting: 


Rp = Rt.sp^. (5) 

The Kepler planetary candidates para meters are esti¬ 
mated assuming the sty is a single star (iBatalha et al.l 
l20III : lBurke et al.ll20l4lMullallv et al.ll20I4ll . and, there- 

fore, may incorrectly report the planet radius if the stel¬ 
lar host is really a multiple star system. The extra flux 
contributed by the companion stars will dilute the ob¬ 
served transit depth, and the derived planet radius de¬ 
pends on the size of star presumed to be transited. The 
ratio of the true planet radius, Rp{true), to the observed 
planet radius assuming a single star with no companions, 
Rp{observed), can be described as: 




Rp{true) 

Rp{observed) 



( 6 ) 


where i?i* is the radius of the (assumed single) primary 
star, and Ft and Rt* are the brightness and the radius, 
respectively, of the star being transited by the planet. 

This ratio reduces to unity in the case of a single star 
(i?t* = i?i* and Ftotai = Ft). For a multiple star system 
where the planet orbits the primary star (i?t* = i?i*), 
the planet size is underestimated only by the flux dilution 
factor: 


Xr = 


Rp {true) 
Rp{observed) 



(7) 


However, if the planet orbits one of the companion stars 
and not the primary star, then the ratio of the primary 
star radius (i?i*) to the radius of the companion star be¬ 
ing transited {Rt*) affects the observed planetary radius, 
in addition to the flux dilution factor. 


3. POSSIBLE COMPANIONS FROM ISOCHRONES 

To explore the possible effects of the undetected stel¬ 
lar companions on the derived planetary parameters, we 
first assess what companions are possible for each KOI. 
For this work, we have downloaded the Cumulative Ke¬ 
pler candidate list and stellar parameters table from the 
NASA Exoplanet Archive. The cumulative list is up¬ 
dated with each new release of the KOI list^; as a result, 
the details of any one star and planet may have changed 
since the analysis for this paper was done. However, the 

® The 2014 October 23 update to the cumula¬ 
tive KOI table was used in the work presented here: 
http: //exoplanetarchive.ipac.caltech.edu/ 
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Figure 1. Radii corrections factors (Xu) are plotted as a function 
of companion-to-primary brightness ratios {bottom axis) and mass 
ratios {top axis) for possible binary systems {top plot) or triple 
{bottom plot) systems. This figure is an example for the G-dwarf 
KOI-299; similar calculations have been made for every KOI. In 
each plot, the dark blue stars represent the correction factors if 
the planet orbits the primary star (equation [TJ; the red circles 
represent the correction factors if the planet orbits the secondary 
star, and the light blue triangles represent the correction factors if 
the planet orbits the tertiary star (equation [6j. The lines are third 
order polynomials fit to the distributions. Unity is marked with a 
horizontal dashed line. 

overall results of the paper presented here should remain 
largely unchanged. 

For the KOI lists, the stellar parameters for each 
KOI were determined by fitting photometric colors and 
spectroscopically derived parameters (wher e available) to 
the Dartmouth Stellar E volution Database (|Dotter et al.l 
120081 iHuber et al.l I20I4I1 . The planet parameters were 
then derived based upon the transit curve htting and 
the associated stellar parameters. Other stars listed in 
the Kepler Input Catalog or UKIRT imaging that may 
be blended with the KOI host stars were accounted for in 
the transit fitting, but, in general, as mentioned above, 
each planetary host star was assumed to be a single star. 

We have restricted the range of possible bound stellar 
companions to each KOI host star by utilizing the same 
Dartmouth isochrones used to determine the stellar pa¬ 
rameters. Possible gravitationally bound companions are 
assumed to lie along on the same isochrone as the pri¬ 
mary star. For each KOI host star, we found the single 
best fit isochrone (characterized by mass, metallicity, and 
age) by minimizing the chi-square ht to the stellar pa¬ 
rameters (effective temperature, surface gravity, radius, 
and metallicity) listed in the KOI table. 


Figure 2. This figure is the same as Figure[2 but for the M-dwarf 
KOI-1085, demonstrating that the details of the derived correction 
factors are dependent upon the KOI properties. 

We did not try to re-derive stellar parameters or in¬ 
dependently find the best isochrone fit for the star; we 
simply identified the appropriate Dartmouth isochrone 
as used in the det ermination of the stellar parameters 
(|Huber et al.ll20ill . We note that there exists an ad¬ 
ditional uncertainty based upon the isochrone finding. 
In this work, we did not try to re-derive the stellar pa¬ 
rameters of the host stars, but rather, we simply hnd 
the appropriate isochrone that matches the KOI stellar 
parameters. Thus, any errors in the stellar parameters 
derivations in the KOI list are propagated here. This is 
likely only a significant source of uncertainty for nearly 
equal brightness companions. 

Once an isochrone was identified for a given star, all 
stars along an isochrone with (absolute) Kepler magni¬ 
tudes fainter than the (absolute) Kepler magnitude of the 
host star were considered to be viable companions; i.e., 
the primary host star was assumed to be the brightest 
star in the system. The fainter companions listed within 
that particular isochrone were then used to establish the 
range of possible planetary radii corrections (equation 
[6|) assuming the host star is actually a binary or triple 
star. Higher order (e.g., quadruple) stellar multiples are 
not considered her e as they represent on ly ^ 3% of the 
stellar population (iRaghavan et ahlDOlOD . 

We have considered six specihc multiplicity scenarios: 

1. single star = 1.0) 

2. binary star - planet orbits primary star 

























4 




8000 7000 6000 5000 4000 3000 8000 7000 6000 5000 4000 3000 

Stellar [K] Stellar [K] 




Figure 3. Top: The mean correction factor (X^) for each KOI is displayed as a function of the effective temperature of the primary star 
(see section 0. The black curves are 3’'^-order polynomials fitted to the distributions (equation [Sj. Bottom: Histograms of the correction 
factors displayed in the top panels. The vertical dashed lines mark the medians of the distributions. Left: The correction factors are 
computed for the Kepler Cumulative Kepler Objects of Interest list. Right: The corrections are computed for the KOIs but assuming 
the KOIs are 10 times closer as may be the case for K2 and TESS. The red points and histograms assume each KOI is single as is the 
case for the published KOI list; the blue points and histograms assume that each KOI has been vetted with radial velocity (RV) and high 
resolution imaging (HRI), and all stellar companions with orbital periods of 2 years or shorter and all stellar com panions located at angular 
distances of > 0.1'' have been detected and accounted for in the planetary radii determinations (see Section 14:^ . For the vetted stars, the 
correction factors are only for undetected stellar companions; detected companions have been assumed to be accounted for in the planet 
radii determinations. 


3. binary star - planet orbits secondary star 

4. triple star - planet orbits primary star 

5. triple star - planet orbits secondary star 

6. triple star - planet orbits tertiary star. 

Based upon the brightness and size differences between 
the primary star and the putative secondary or tertiary 
companions, we have calculated for each KOI the pos¬ 
sible factor by which the planetary radii are underesti¬ 
mated (Xji). If the star is single, the correction factor is 
unity, and if, in a multiple star system, the planet orbits 
the primary star, only flux dilution affects the observed 
transit depth and the derived planetary radius (eq. [7]). 

For the scenarios where the planets orbits the sec¬ 
ondary or tertiary star, the planet size correction fac¬ 
tors (eq. [6]) were determined only for stellar companions 
where the stellar companion could physically account for 


the observed transit depth. If more than 100% of the stel¬ 
lar companion light had to be eclipsed in order to produce 
the observed transit in the presence of the flux dilution, 
then that star (and all subsequent stars on the isochrone 
with lower mass) was not considered viable as a poten¬ 
tial source of the transit. For example, for an observed 
1% transit, no binary companions can be fainter than 
the primary star by 5 magnitudes or more; an eclipse of 
such a secondary star would need to be more than 100% 
deep. The stellar brightness limits were calculated inde¬ 
pendently for each planet within a KOI system so as to 
not assume that all planets within a system necessarily 
orbited the same star. 

Figures [Hand [7] show representative correction factors 
(Xji) for KOI-299 (a G-dwarf with a super-Earth sized 
Rp = 1.8 ± O.24ii0 planet) and for KOT1085 (an M- 
dwarf with an Earth-sized Rp = 0.92 ± O.lSi?® planet). 
The planet radius correction factors {Xu) are shown as 
a function of the companion-to-primary brightness ratio 
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(bottom x-axis of plots) and the companion-to-primary 
mass ratio (top x-axis of plots) and are determined for 
the KOI assuming it is a binary-star system (top plot) 
or a triple-star system (bottom plot). 

The amplitude of the correction factor (X/j) varies 
strongly depending on the particular system and which 
star the planet may orbit. If the planet orbits the pri¬ 
mary star, then the largest the correction factors are for 
equal brightness companions {^/2 ~ 1.4 for a binary sys¬ 
tem and -x/S ^ 1.7 for a triple system) with an asymp¬ 
totic approach to unity as the companion stars become 
fainter and fainter. If the planet orbits the secondary or 
tertiary star, the planet radius correction factor can be 
significantly larger - ranging from ^2 — 5 for binary 
systems and X/j > 2 — 20 or more for triple systems - 
depending on the size and brightness of the secondary or 
tertiary star. 

4. MEAN RADII CORRECTION FACTORS {Xr) 

It is important to recognize the full range of the possi¬ 
ble correction factors, but in order to have a better under¬ 
standing of the statistical correction any given KOI (or 
the KOI list as a whole) may need, we must understand 
the mean correction for any one multiplicity scenario and 
convert these into a single mean correction factor for each 
star. To do this, we must take into account the proba¬ 
bility the star may be a multiple star, the distribution 
of mass ratios if the star is a multiple, the probability 
that the planet orbits any one star if the stellar system 
has multiple stars, and whether or not the star has been 
vetted (and how well it is has been vetted) for stellar 
companions. 

In order to calculate an average correction factor for 
each multiplicity scenario, we have fitted the individual 
scenario correction factors as a function of mass ratio 
with a 3’''^-order polynomial (see Fig. [T]and[2|). Because 
the isochrones are not evenly sampled in mass, taking a 
mean straight from the isochrone points would skew the 
results; the polynomial parameterization of the correc¬ 
tion factor as a function of the mass ratio enables a more 
robust determination of the mean correction factor for 
each multiplicity scenario. 

If the companion-to-primary mass ratio distribution 
was uniform across all mass ratios, then a straight mean 
of the correction values determined from each polynomial 
curve would yield the average correction factor for each 
multiplicity scenario. However, the mass ratio distribu¬ 
tion is likely not uniform, and we have adopted t he form 
displayed in Figure 16 of iRaehavan et all (120101) . That 
distribution is a nearly-flat frequency distribution across 
all mass ratios with a ~ 2.5x enhancement for nearly 
equal mass companion stars {q > 0.95). This distribu- 
tion is in contrast to t he Ga ussian distribution shown in 
iDuouenn ov fc Mavoil (II 99Ill: ho wever, the more recent 
results of iRaghavan et al.l (|2010ll incorporate more stars, 
a broader breadth of stellar properties, and multiple com¬ 
panion detection techniques. 

The mass ratio distribution is convolved with the poly¬ 
nomial curves fitted for each multiplicity scenario, and a 
weighted mean for each multiplicity scenario was calcu¬ 
lated for every KOI. For example, in the case of KOI-299 
(Fig. [T|), the single star mean correction factor is I.O (by 
definition). For the binary star cases, the average sce¬ 
nario correction factors are 1.14 (planet orbits primary) 


and 2.28 (planet orbits secondary); for the triple stars 
cases, the correction factors are 1.16 (planet orbits pri¬ 
mary), 2.75 (planet orbits secondary), and 4.61 (planet 
orbits tertiary). For KOl-1085 (Fig. [2|), the weighted 
mean correction factors are 1.18, 1.56, 1.24, 1.61, and 
2.29, respectively. 

To turn these individual scenario correction factors 
into an overall single mean correction factor {Xu) per 
KOI, the six scenario corrections are convolved with the 
probability that a KOI will be a single star, a binary star, 
or a triple star. The multiplicity ra te of the Kepler stars 
is still unclear (IWang et al.l 12014^ . and, indeed, there 
may be some contradictory evidence for the the exact 
value for the mult i plicity rates of the K OI host stars (e.g., 
iHorch et al.ll20l4lWang et al.ll20I4bl) . but the multiplic¬ 
ity rates appear to be near 40%, similar to the general 
held population. In the absence of a more dehnitive esti- 
mate , we have chosen to uti lize the multiplicity fractions 
from iRaghavan et al.l (|20IClll : a 54% single star fraction, 
a 34% binary star fract ion, and a 12% triple star fraction 
(|Raghavan et al.l [201011 . We have grouped all higher or¬ 
der multiples (3-I-) into the single category of “triples”, 
given the relatively rarity of the quadruple and higher 
order stellar systems. For the scenarios where there are 
multiple stars in a system, we have assumed that the 
planets are equally likely to orbit any one of the stars 
(50% for binaries, 33.3% for triples). 

The hnal mean correction factors (X^) per KOI are 
displayed in Figure [3] the median value of the correc¬ 
tion factor and the dispersion around that median is 
(X/j) = 1.49 ± 0.12. This median correction factor im¬ 
plies that assuming a star in the KOI list is single, in the 
absence of any (observational) companion vetting, yields 
a statistical bias on the derived planetary radii where the 
radii are underestimated, on average, by a factor ^ 1.5, 
and the mass density of the planets are overestimated by 
a factor of ~ 1.5^ ^ 3. 

From Figure [31 it is clear that the mean correction 
factor {Xu) depends upon the stellar temperature of the 
host star. As most of the stars in the KOI list are dwarfs, 
the lower temperature stars are typically lower mass stars 
and, thus, have a smaller range of possible stellar com¬ 
panions. Thus, an average value for the correction factor 
1.5 represents the sample as a whole, but a more accurate 
value for the correction factor can be derived for a given 
star, with a temperature between 3500 7500K, 

using the htted 3''‘^-order polynomial: 

{Xr) = a3{Teff)^ + a 2 {Teff)‘^ + ai{Teff) + ao (8) 

where 03 = —1.19118 x 10“^^, 02 = 1.61749 x 10“^, oi = 
—0.000560, and oq = 1.64668. 

In the absence of any specific knowledge of the stel¬ 
lar properties (other than the effective temperature) and 
in the absence of any radial velocity or high resolution 
imaging to assess the specific companion properties of a 
given KOI, (see section 02]), the above parameterization 
(equation |5|) can be used to derive a mean radii correc¬ 
tion factor (Xfl) for a given star. For G-dwarfs and hot¬ 
ter stars, the correction factor is near (X/j) ~ 1.6. As 
the stellar temperature (mass) of the primary decreases 
to the range of M-dwarfs, the correction factor can be as 
low as {Xji) ^ 1 . 2 . 
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Figure 4. The distribution of the ratio of the total planetary radii 
uncertainties (crpi ^) to the quoted radii uncertainties )from the 

cumulative KOI list (see equation [Toll. For the red histogram, it is 
assumed that the KOIs are single as is the case in the published 
KOI list; for the blue histogram, it is assumed that each KOI has 
been vetted with radial velocity (RV) and high resolution imaging 
(see Section (Mil. The vertical dashed lines represent the median 
values of the distributions: {’^Rp/^'up) ~ the unvetted 

KOIs and { oRplo'ji ) = 1-15 for the vetted KOIs (see section [4. 2 II . 

4.1. Planet Radius Uncertainty Term from {Xu) 

The mean correction factor is useful for understand¬ 
ing how strongly the planetary radii may be underesti¬ 
mated, but an additional uncertainty term derived from 
the mean radius correction factor is potentially more 
useful as it can be added in quadrature to the formal 
planetary radii uncertainties. The formal uncertainties, 
presented in the KOI list, are derived from the uncer¬ 
tainties in the transit fitting and the uncertainty in the 
knowledge of the stellar radius, and they are calculated 
assuming the KOIs are single stars. We can estimate 
an additional planet radius uncertainty term based upon 
the mean radii correction factor as 

= \{Xr)Rp - Rp\ = \{Xr) - 1.0\Rp (9) 

where Rp is the observed radius of the planet. Adding in 
quadrature to the reported uncertainty, a more complete 
uncertainty on the planetary radius can be reported as 

/ \ 1/2 

= + ( 10 ) 

where a'j^^ is the uncertainty of the planetary radius as 
presented in the KOI list. 

The distribution of the ratio of the more complete KOI 
radius uncertainties {o'Rp) to the reported KOI radius 
uncertainties (tr^j ) is shown in Figured Including the 
possibility that a KOI may be a multiple star increases 
the planetary radii uncertainties. While the distribution 
has a long tail dependent upon the specific system, the 
planetary radii uncertainties are underestimated as re¬ 
ported in the KOI list, on average, by a factor of 1.7. 

4.2. Effectiveness of Companion Vetting 

The above analysis has assumed that the KOIs have 
undergone no companion vetting, as is the assumption 
in the current KOI list. In reality, the Kepler Project 


has funded a substantial ground-based follow-up obser¬ 
vation program which includes radial velocity vetting and 
high resolution imaging. In this section, we explore the 
effectiveness of the observational vetting. 

The observational vetting reduces the fraction of unde¬ 
tected companions. If there is no vetting or all stars are 
assumed to be single, as is the case for the published KOI 
list, then the fraction of undetected companions is 100% 
and the mean correction factors {Xr) are as presented 
above. If every stellar companion is detected and ac¬ 
counted for in the planetary parameter derivations, then 
the fraction of undetected companions is 0%, and the 
mean correction factors are unity. Reality is somewhere 
in between these two extremes. 

To explore the effectiveness of the observational vetting 
on reducing the radii corrections factors (and the asso¬ 
ciated radii uncertainties), we have assumed that every 
KOI has been vetted equally, and all companions within 
the reach of the observations have been detected and ac¬ 
counted. Thus, the corrections factors depend only on 
the fraction of companions stars that remain out of the 
reach of vetting and undetected. 

In this simulation, we have assumed that all compan¬ 
ions with orbital periods of 2 years or less and all compan¬ 
ions with angular separations of 0.1" or greater have been 
detected. This, of course, will not quite be true as ran¬ 
dom orbital phase effects, inclination effects, companion 
mass distribution, stellar rotation effects, etc. will dimin¬ 
ish the efficiency of the observations to detect compan¬ 
ions. We recognize the simplicity of these assumptions; 
however, the purpose of this section is to assess the use¬ 
fulness of observational vetting on reducing the uncer¬ 
tainties of the planetary radii estimates, not to explore 
fully the sensitivities and completeness of the vetting. 

Typical follow-up observations include stellar spec¬ 
troscopy, a few radial velocity measurements, and high 
resolution imaging. The radial velocity observations usu¬ 
ally include 2 — 3 measurements over the span of 6 — 9 
months and are typically sufficient to identify potential 
stellar companions with orbital periods of ;< 1 — 2 years or 
less. While determining full orbits and stellar masses for 
any stellar companions detected typically requires more 
intensive observing, we have estimated that 3 measure¬ 
ments spanning 6 — 9 months is sufficient to enable the 
detection of an RV trend for orbital periods of ~ 2 years 
or less and mark the star as needing more detailed obser¬ 
vations. The amplitude of the RV signature, and hence 
the ability to detect companions, does depend upon the 
masses of the primary and companion stars; massive 
stars with low mass companions will display relatively 
low RV signatures. However, RV vetting for the Kepler 
program has been done at a level of < 100 — 200 m/s, 
which is sufficient to detect (at > 4 — 5cr) a late-type 
M-dwarf companion in a two-year orbit around a mid 
B-dwarf primary. Indeed, the RV vetting is made even 
more effect ive by searching f or companions via spectral 
signatures (iKolbl et al.ll2015il . 

The high resolution imaging via adaptive optics. 


has resolutions of 0.02" — 0.1" (e.g.. iHowell et al.ll201ll: 

Horch et al. 

120121 ILillo-Box et al.l 120121: [Adams et al. 

20121 2013 

: Dressing et al.1 2014 [Law et al.l 2014: 

Lillo-Box et al.l 2014: Wane et al. 2014a: Everett et al. 

2014 iHorch 

et al.ll20l4 [Gilliland et al.ll20l4 [Star et al. 
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Figure 5. Top Left: Example of the stellar companion period distribution that is vetted by radial velocity (RV) monitoring and high 
resolution i maging (HRI) for KOI- 299. The Gaussian curve, normalized to unity, is the log-normal orbital period distribution of stellar 
companions IIRaghavan et ani20inii , and the hatched regions mark where the potential observational vetting is assumed to have detected all 
companions in that period range. The solid (red) region of the Gaussian cor responds to the fraction of companions (31% for KOT299) that 
would remain undetected by the RV and HRI observations (see section EM- Bottom Left: The distribution of the fraction of companions 
across all KOIs that remain undetected by observational radial velocity and high resolution imaging vetting (i.e., the red area in the top 
figure but for all KOIs). The vertical dashed line represents the median fraction (41%) of companions left undetected by the observational 
vetting of all KOIs. Top Right: The same as the top left figure but assuming that KOI-299 is 10 times closer as may be the case for 
K2/TESS targets. At such a close distance for a hypothetical KOI-299, RV and HRI vetting leaves only about 7% of the possible stellar 
companions undetected. Bottom Right: As the bottom left figure but assuming all of the KOI host stars are 10 times closer as may be the 
case for K2/TESS. The vertical dashed line represents the median fraction (16%) of companions left undetected if the KOIs were at these 
closer distances. 


I2014D . Based upon Monte Carlo simulations in which 
we have averaged over random orbital inclinations and 
eccentricities, we have calculated the fraction of time 
within its orbit a companion will be detectable via high 
resolution imaging. With typical high resolution imag¬ 
ing of 0.05", we have estimated that > 50% of the stel¬ 
lar companions will be detected at one full-width half¬ 
maximum (FWHM=0.05") of the image resolution and 
beyond and > 90% at > 2 FWHM (0.1") of the image 
resolution and beyond. 

To determine what fraction of possible stellar compan¬ 
ions would be detected in such a scenario, we have used 
the nearly log-normal orbital period distribution from 
iRaghavan et al.l (|2010H . To convert the high resolution 
imaging limits into period-limits, we have estimated the 
distance to each KOI by determining a distance modulus 
from the observed Kepler magnitude and the absolute 


Kepler magnitude associated with the fitted isochrone. 
The median distance to the KOIs was found to be ^ 900 
pc, corresponding to ^ 90 AU for O.I" imaging. Using 
the isochrone stellar mass, the semi-major axis detection 
limits were converted to orbital period limits (assuming 
circular orbits). 

Combining the 2-year radial velocity limit and the O.I" 
imaging limit, we were able to estimate the fraction of 
undetected companions for each individual KOI (see Fig¬ 
ure 0. The distribution of the fraction of undetected 
companions ranges from ~ 20 — 60% and, on average, 
the ground-based observations leave « 41% of the possi¬ 
ble companions undetected for the KOIs (see Figure 0. 

The mean correction factors (Xu) are only applicable 
to the undetected companions. For the stars that are 
vetted with radial velocity and/or high resolution imag¬ 
ing, the intrinsic stellar companion rate for the KOIs of 































46% (|Raghavan et aljr2010r) is reduced by the unvetted 
companion fraction for each KOI. That is, we assume 
that companion stars detected in the vetting have been 
accounted for in the planetary radii determinations, and 
the unvetted companion fraction is the relevant compan¬ 
ion rate for determining the correction factors. In the 
KOI-299 example (Fig. [5]), the undetected companion 
rate used to calculate the mean radii correction factor is 
0.46 X 0.31 = 0.1426. This lower fraction of undetected 
companions in turn reduces the mean correction factors 
for the vetted stars which are displayed in Figure [3] (blue 
points). 

Instead of a mean correction factor of (Xu) ^1.5, the 
average correction factor is (Xu) = 1.20 ± 0.06 if the 
stars are vetted with radial velocity and high resolution 
imaging. The mean correction factor still changes as a 
function of the primary star effective temperature but 
the dependence is much more shallow with coefficients 
for equation[S]of 03 = —6.73847 x 10“^^, 02 = 9.38966 x 
10“®,ai = —0.000352, and uq = 1.40391 (see Figure [S]). 

4.3. K2 and TESS 

The above analysis has concentrated on the Kepler 
Mission and the associated KOI list, but the same 
effects will apply to all transit surveys inc l uding K2 
(|Howell et al.ll20lll and TESS (iRicker et al.ll20lj[i . If 
the planetary host stars from K2 and TESS are also 
assumed to be single with no observational vetting, 
the planetary radii will be underestimated by the same 
amount as the Kepler KOIs (Fig. [3] and eq. [ 8 ]). 

Many K2 targets and nearly all of TESS targets will be 
stars that are typically 4—5 magnitudes brighter than the 
stars observed by Kepler, and therefore, K2 and TESS 
targets will be ^ 10 times closer than the Kepler targets. 
The effectiveness of the radial velocity vetting will re¬ 
main mostly unaffected by the brighter and closer stars, 
but the effectiveness of the high resolution imaging will 
be significantly enhanced. Instead of probing the stars 
to within ^ lOOAU, the imaging will be able to detect 
companion stars within ^ 10 AU of the stars. 

As a result, the fraction of undetected companions will 
decrease significantly. Even for the Kepler stars that 
undergo vetting via radial velocity and high resolution 
imaging, ~ 40% of the companions remain undetected. 
But for the stars that are 10 times closer that fraction 
decreases to ~ 16% (see Figure (H). This has the strong 
benefit of greatly reducing the mean correction factors 
for the stars that are observed by K2 and TESS and 
are vetted for companions with radial velocity and high 
resolution imaging. 

The mean correction factor for vetted K2/TESS-like 
stars is only {Xu) = 1.07 ± 0.03. The correction fac¬ 
tor has a much flatter dependence on the primary star 
effective temperature, because the majority of the possi¬ 
ble stellar companions are detected by the vetting. The 
coefficients for equation [S] become 03 = —4.12309 x 
10-12^02 = 5.89709 X 10"®, m = -0.000242, and oq = 
1.30060. The mean radii correction factors for vetted 
K2/TESS planetary host stars correspond to a correc¬ 
tion to the planetary radii uncertainties of only ~ 2 %, 
in comparison to a correction of ~ 70% if the K2/TESS 
stars remain unvetted. 

For K2 and TESS, where the number of candidate 
planetary systems may outnumber the KOIs by an or¬ 


der of magnitude (or more), single epoch high resolution 
imaging may prove to be the most important observa¬ 
tional vetting performed. While the imaging will not 
reach the innermost stellar companions, radial velocity 
observations require multiple visits over a baseline com¬ 
parable to the orbital periods an observer is trying to 
sample. In contrast, the high resolution imaging requires 
a single visit (or perhaps one per filter on a single night) 
and will sample the majority of the expected stellar com¬ 
panion period distribution. 

5. EFFECT OF UNDETECTED COMPANIONS ON THE 
DERIVED OCCURRENCE RATES 

Understanding the occurrence rates of the Earth-sized 
planets is one of the primary goals of the Kepler missio n 
and one of the uses of the KOI list (jBorucki et al.l[ 20 T(ill . 
It has been shown that the transition from ’’rocky” to 
”non-rocky” planets occurs near a radius of R „ = 1.6 
and the transition is very sharp (|Rogersll2014ll . However, 
the amplitude of the uncertainties resulting from unde¬ 
tected companions may be large enough to push plan¬ 
ets across this boundary and affect our knowledge of the 
fraction of Earth-sized planets. 

We have explored the possible effects of undetected 
companions on the derived occurrence rates. The plan¬ 
etary radii can not simply be multiplied by a mean cor¬ 
rection factor Xji, as that factor is only a measure of the 
statistical uncertainty of the planetary radius resulting 
from assuming the stars are single and only a fraction 
of the stars are truly multiples. Instead a Monte Carlo 
simulation has been performed to assign randomly the 
effect of unseen companions on the KOIs. 

The simulation was performed 10,000 times for each 
KOI. For each realization of the simulation, we have 
randomly assigned the star to be single, binary, or 
triple star via the 54% . the 34% and the 12% fractions 
(|Raghavan et al.ir2010D . If the KOI is assigned to be a 
single star, the mean correction factors for the planets 
in that system are unity: (Xji) = 1. If the KOI star 
is a multiple star system, we have randomly assigned 
the stellar companion masses according to the masses 
available from the fi tted isochrones and us ing the mass 
ratio distribution of (|Raghavan et al.l[2010H . Finally, the 
planets are randomly assigned to the primary or to the 
companion stars (i.e., 50% fractions for binary stars and 
33.3% fractions for triple stars). Once the details for the 
system are set for a particular realization, the final cor¬ 
rection factor for the planets are determined from the 
polynomial fits for the individual multiplicity scenarios 
(e.g.. Fig. [Hand[2]). 

For each set of the simulations, we compiled the frac¬ 
tion of planets within the following planet-radii bins: 
Rp < 1.6 1.6 < Rp < 3.0 Rp] 3.0 < i?p < 10 R^ cor¬ 

responding to Earth-sized, super-Earth/mini-Neptune- 
sized, and Neptune-to-Jupiter-sized planets. The raw 
fractions directly from the KOTlist, for these three cat¬ 
egories of planets, are 33.3%, 46.0%, and 20.7%. Note 
that these are the raw fractions and are not corrected for 
completeness or detectability as must be done for a true 
occurrence rate calculation; these fractions are necessary 
for comparing how unseen companions affect the deter¬ 
mination of fractions. Finally, we repeated the simula¬ 
tions, but using the undetected multiple star fractions af¬ 
ter vetting with radial velocity and high resolution imag- 
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Figure 6. The percent change of the measured occurrence rates from the raw KOI list (i.e., uncorrected for completeness) caused by 
undetected stellar companions and assuming all of the stars in the published KOI list are single. The fractional changes are computed 
for three different planet size bins as labeled in the plots. The results of assuming all the stars are single (i.e ., un vetted) are in the red 
histogram and the results of vetting the all of the stars with radial velocity and high resolution imaging (section l4.2ll are shown in the blue 
histogram. The histograms are based upon Monte Carlo simulations described in section [5] 


ing had been performed, thus, effectively increasing the 
fraction of stars with correction factors of unity. 

The distributions of the change in the fractions of plan¬ 
ets in each planet category, compared to the raw KOI 
fractions, are shown in Figure |6l If the occurrence rates 
utilize the assumed-single KOI list (i.e., unvetted), then 
the Earth-sized planet fraction may be overestimated by 
as much as 15 — 20% and the giant-planet fraction may 
be underestimated by as much as 30%. Interestingly, the 
fraction of super-Earth/mini-Neptune planets does not 
change substantially; this is a result of smaller planets 
moving into this bin, and larger planets moving out of 
the bin. In contrast, if all of the KOIs undergo vetting 
via radial velocity and high resolution imaging, the frac¬ 
tional changes to these bin fractions are much smaller: 
5 — 7 % for the Earth-sized planets and 10 — 12% for the 
Neptune/Jupiter-sizes planets. 

6. SUMMARY 

We present an exploration of the effect of undetected 
companions on the measured radii of planets in the Ke¬ 
pler sample. We find that if stars are assumed to be 
single (as they are in the current Kepler Objects of In¬ 
terest list) and no companion vetting with radial velocity 
and/or high resolution imaging is performed, the plane¬ 
tary radii are underestimated, on average, by a factor of 
(Xji) = 1.5, corresponding to an overestimation of the 
planet bulk density by a factor of ~ 3. Because lower 
mass stars will have a smaller range of stellar companion 
masses than higher mass stars, the planet radius mean 
correction factor has been quantified as a function of stel¬ 
lar effective temperature. 

If the KOIs are vetted with radial velocity observations 
and high resolution imaging, the planetary radius mean 
correction necessary to account for undetected compan¬ 
ions is reduced significantly to a factor of (Xu) = 1 . 2 . 
The benefit of radial velocity and imaging vetting is even 
more powerful for missions like K2 and TESS, where the 
targets are, on average, ten times closer than the Kepler 


Objects of Interest. With vetting, the planetary radii for 
K2 and TESS targets will only be underestimated, on 
average, by 10%. Given the large number of candidates 
expected to be produced by K2 and TESS, single epoch 
high resolution imaging may be the most effective and 
efficient means of reducing the mean planetary radius 
correction factor. 

Finally, we explored the effects of undetected compan¬ 
ions on the occurrence rate calculations for Earth-sized, 
super-Earth/mini-Neptune-sized, and Neptune-sized and 
larger planets. We find that if the Kepler Objects of In¬ 
terest are all assumed to be single (as they currently are 
in the KOI list), then the fraction of Earth-sized plan¬ 
ets may be overestimated by as much as 15-20% and the 
fraction of large planets may be underestimated by as 
much as 30% 

The particular radial velocity observations or high res¬ 
olution imaging vetting that any one KOI may (or may 
not) have undergone differs from star to star. Companion 
vetting simulations presented here show that a full un¬ 
derstanding and characterization of the planetary com¬ 
panions is dependent upon also understanding the pres¬ 
ence of stellar companions, but is also dependent upon 
understanding the limits of those observations. For a 
final occurrence rate determination of Earth-sized plan¬ 
ets and, more importantly, an uncertainty on that oc¬ 
currence rate, the stellar companion detections (or lack 
thereof) must be taken into account. 

The authors would like to thank Ji Wang, Tim Mor¬ 
ton, and Gerard van Belle for useful discussions during 
the writing of this paper. This research has made use of 
the NASA Exoplanet Archive, which is operated by the 
California Institute of Technology, under contract with 
the National Aeronautics and Space Administration un¬ 
der the Exoplanet Exploration Program. Portions of this 
work were performed at the California Institute of Tech¬ 
nology under contract with the National Aeronautics and 
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